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Abstract 

We derive the asymptotical dynamical law for Ginzburg-Landau vortices in an inhomo- 
geneous background density under the Schrodinger dynamics, when the Ginzburg-Landau 
parameter goes to zero. New ingredients involve across the cores lower bound estimates 
and approximations. 

1 Introduction 

We are interested in the two dimensional Gross-Pitaevskii equation 

idtu- Au + ^{V{x) + \u\^)u = (GP) 

for n : X R+ C, where < e ^ 1 and 1/ : M — t- ]R+ is a smooth potential such that 

V{x) — -|-oo as |x| — )■ -|-oo. 

The Gross-Pitaevskii equation is a widely used model to describe the dynamics of a Bose- 
Einstein condensate in a trapping potential V. The equation on arises via dimension 
reduction from 3 dimensions; this has been justified for particular choices of y in [1] for 
example. 

Equation ()GPp is Hamiltonian, with Hamiltonian given by 

VnP 1 , luP lul'' 



Another quantity which is preserved by the flow associated with (1GP|) is the total mass M, 
given by 

M(n) = / |m|2. 



For each m > 0, there exist£] at least one positive ground state t] = rji.^m '■ — ^ I^"*" of total 
mass equal to m. By definition, a ground state 77 realizes the infimum 

£e,v{v) = mf{£,y{g), g G H\r\C), M{g) = m}, 



We refer to Section [9] for the details on a number of statements regarding the ground states which we 
state without justification in this introduction. 
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and satisfies the Euler-Lagrange equation 

1 1 

-Ari + ^{V + ri^)r] = —At?, 

where we write the Lagrange multipHer as for some A = A^^m- 
In the hmit e — t- 0, we have 

rf ^ Ptf inL2(M2)^ 

where the function p^F, known as the Thomas-Fermi profile in the physics hterature, is given 
by Ptf{x) := {Xo — V)^{x) where the number Aq is uniquely determined by the mass condition 

/ {Xq — V{x))'^ dx = m. 

We will study the behaviour of solutions of (jGPp which correspond, in a sense to be made 
precise later, to perturbations of the ground state r? by a finite number of quantized vortices, 
each carrying a single quantum of vorticity. Our goal is to prove that these vortices persist, 
and to describe their evolution in time. 

We will show that to leading order the vortices do not interact, and that each one evolves 
(in a renormalized time scale) by the orthogonal gradient flow for the function log Ptf i with 
a sign depending on the winding number of the given vortex. More precisely, let 

^TF '■= {X ■ Ptf{x) > 0} 

be the interior of the limiting supporlH of the ground state, let {b^}l^i be distinct points in 
fi-TF) and let di, . . . ,di G {—1, +!}• For each i G {1, • • • , /}, we denote by the solution of 
the ordinary differential equation 

bS) = d,^^{b,{t)), (1) 

Ptf 

where V"*- = {—dx2,dx^), with initial datum bi{0) = b^. 

Theorem 1. Let (ng)e>o be a family of initial data for (|GPp such that 

M(u°) = m, 
I 

£s,v{u^e) < £e,viri)+T^^PTFib'',)\loge\ + o(|loge|), 

i=l 

and ^ 

curl(^) 27r V (i,5,o in W-^'^Utf), 
Ptf .^^ ^ 

as e — 7- 0. Then, as long as the points {bi(t)}[^i remain distinct, 

t\ ' 

curl(^) ^ 27r di6h^(t) in W~J'\nTF), 
Ptf .^^ 

as e — )• 0, where := Us{-,t\loge\). 

^Note that we do not assume that Qtf is simply connected or that its boundary is smooth. 
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Here, := (iu*, Vn*) where {z,w) := lm{zw). Therefore, 



curl(^) = -curl J (-^) = 

is the Jacobian determinant of u\l ^ p^F- It is widely recognized, in the present regime for 
the Ginzburg-Landau energy, that the notion of a vortex of winding number di located at 
the point 6j(t) is appropriately described by the presence of the term 2TTdi6i,.(^i-^ in the limit 
of the vorticity field curl j[ul/y/p^). 

Remark 1. Note that the ordinary differential equations ([1]) are decoupled. Also, since 
PTF{bi{t)) = pTF{b^) for any t G M the points {bi{t)}\^i remain distinct for all times unless 
two of them are located on the same level line of Ptf and have opposite circulations. 

Results of this sort in the homogeneous case r/ = 1 were first proved in the late 1990s, 
see [5l [6], and have subsequently been developed by a number of authors, see for example 
[151 [21 [Tl]. The point of this paper is thus to understand the effect of the inhomogeneity on 
the dynamical law for the vortices. 

We remark that a number of authors have studied questions about vortex dynamics in in- 
homogeneous backgrounds for parabolic equations [14:\ [T2] , or more recently [19] for a quite 
general class of equations of mixed parabolic-Schrodinger type. The case of pure Schrodinger 
dynamics presents distinct difficulties and as far as we know has not been treated until now. 

The sequel of this introduction is devoted to the presentation of the strategy leading to 
Theorem [1] We notice that will actually prove a result (Theorem [2] below) which is stronger 
in two respects than Theorem [1) first it describes the dynamics of vortices at small but fixed 
value of e, rather than asymptotically as e — )• in Theorem [H and second it applies to a 
broader class of inhomogeneous equations (see (jNHGp below) where rj need not necessarily 
be the profile of a ground state. 

1.1 Perturbation equation and Theorem [2] 

For the class of initial data which we consider in Theorem [H it is convenient to rewrite the 
corresponding solutions of (IGPp in the form 

u{x,t) = r]{x)w{x,t) (2) 

and to study the evolution equation for w. One easily 
then w solves 

irfdiw — div(r/^V'w) H — 2^^(I^P 
In particular, the change of phase and time scale 

leads to the equation 

i|log e\ri^dtv - dlw{ri^Vv) + \r]^{\v\^ - l)v = ^ (NHG) 



checks that if u is a solution to (|GPp , 

— \)W = ttT] W. 
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for V. Note that the change of time scale is related to the fact that the phenomenon which we 
wish to describe, namely vortex motion, arises in times of order one in that new time scale 
(see the definition of u\ in the statement of Theorem [T]) . 

Our analysis will henceforth focus on equation (jNHGp . Equation (jNHGp . like ()GPp . is 
Hamiltonian, with Hamiltonian given by the weighted Ginzburg-Landau energy 

EeM^I esAv)=[ ^^^^ + ^^^^4^- (3) 



2 ' ' 4e2 

As a matter of fact, using the Euler-Lagrange equation for rj, one realizes that 

£e,viu) = £s,viv) + Ee,^{v) + ^ (M(n) - M(r?)) . (4) 

In the sequel, we enlarge our framework and consider equation (jNHGp where t] : — )• M 
is any smooth positive function such that the corresponding Cauchy problem is globally 
well-posed for initial data in {M."^ , r] dx) and such that the corresponding solutions can be 
approximated by smooth solution^. In particular, under those assumptions the energy E^^^^ 
is preserved along the flow of (jNHGp . 

Let e > and let i7 C be a bounded open set. Let {a^}[^i be distinct points in Q, and 
let di, . . . ,di G {—1, +!}• For each i £ {!,■ ■ ■ ,1}, we denote by ai{t) the solution, as long as 
it does not reach dQ, of the ordinary differential equation 

di{t) = diV^ log 7]\ai{t)), (5) 

where V"*- = {—dx2,dx:^), with initial datum aj(0) = a?. 
We assume that 

rjmin ■■= inf r]{x) > 0, 
and we fix a time T^oi > such that 

Pmin ■■= min min {{l-d{ai{t), aj{t))} i^jLl{d{ai{t),dQ)}iLl{l}} > 0. (6) 

Finally, we consider a finite energy solution v of (jNHGp . we set := v{-,t) and we define, 
for t G [0,Teoi], 

I 

1=1 

and 

S*:=%^-vr5:,^(a.(t)). (8) 
|loge| ^ 

We will deduce Theorem [1] from 



^This can be verified for a wide variety of weight functions rj, but we wish not consider that discussion here 
since we already know by means of the change of unknown ([2} that it satisfied when is a ground state. 
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Theorem 2. There exist positive constants eo, 7o and Co, depending only on I, ij^in, Pmin, 
and \\V'ii'^\\l°°[q), such that if < e < Eq and if TP + r^< 70, then 

ri<rl+ (e» + rl + (,« _ j) + (9) 

as long as t < T^ox and TP + r* (t) < 79. 

Remark 2. i) As we shall discuss in Section 11.21 below, the quantity r* , which is a sort 
of discrepancy measure, can be thought of as measuring the distances between the "actual 
vortex locations" and the desired vortex locations. The quantity S*, multiplied by |loge|, 
corresponds to the excess of energy of the solution v with respect to an energy minimizing 
field possessing the vortices at the points ai{t). Notice that since S^^r, is preserved by the flow 
for V and rf is preserved by the flow for the a[s, we have 

= VtG[0,T]. 

a) Theorem [2] is interesting for initial data such that T^ + is small. The existence of 
such data is standard. For example, if we fix / : [0, 00) — )• [0, 1] such that /' > 0, /(O) = 0, 
and /(s) — )• 1 as s — )• 00, then for the initial data 

:= J] /(^^) ( ^~_ \ \ , x:={xi,X2)^z = Xi+ix2, 

one can check that T^ < C|loge|~^, < Ce. In any case, ([9]) contains the error term in 
log |loge|/|loge| which implies that ([9]) only yield the inequality S'^ + r* < 70 for times at 
most of order log |loge|. 

Hi) One could supplement the claims of Theorem [2] with closeness estimate for j{v) to a 
reference field j* of very simple form. This would follow from an application of Corollary [1] 
below; however at the level of approximation which we have adopted here it is only meaningful 
in a neighborhood of size o(l/|loge|) of the vortex core. 

1.2 Elements in the proofs 

Under the conditions that will prevail throughout most of this paper, we will be able to 
identify points ^^i, . . . and a number such that 

I 

\\Jv* -T^Y.di^^\\w-^^^{n) <r\ < e^'^ < r*. (10) 

i=l 

This is expressed in Proposition [2] below, and entitles us to think of i = 1, . . . , Z as being 
the "actual locations of the vortices" in up to precision of order < r^. Admitting this 
interpretation, then basic facts about the norm, recalled in Section [21 imply that 

1 ^ 

r* =-(l + o(l))5]|e*-a*| (11) 
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is essentially the aggregate distance between the actual vortex locations and the desired 
vortex locations, as remarked above. 

Heuristic considerations also suggest that if is a function with vortices at points 
(,1, . . . (or more precisely, if holds), then 

E,,,{v')>7r\loge\{l-o{l))Y,v\Cl). (12) 

i=l 

Hence £'e^,,(f*) — 7r|log e| ^ -^^ ?7^(^*) corresponds to energy that is not committed to the 
vortices, and this energy in principle can cause difficulties for our analysis. From (llOp . (Ilip . 
we have 



,,(t;*)-^|loge|^7?2(e*) < llogel ('s* + i(l + o(l))||V7?2|Ur*] 
i=i V ^ / 



(13) 

< |loge| (So + Cr*). 



Ee 



For our analysis, it suffices to use estimates in the spirit of (jl3p that are a little weaker than 
those suggested in (llSp , these are established in Proposition [TJ We expect from (112p and 
(jl3|) that control of r* should yield a good deal of information about v*. This is expressed in 
Proposition [3l where we compare j{v^) to a reference field j*. An important feature of that 
approximation is that it holds across the vortex core. 

In order to control the evolution in time of r* , we rely on some evolution equations satisfied 
by smooth solutions of (|NHGp . Conservation of energy is a consequence of the identity 

dte.^^iv) = diyir]\Vv,Vt)), (14) 

and the canonical equation for conservation of mass can be written 

^^dt{ri\\v\'-l))=V-{r^'jiv)). (15) 
The vorticity Jv satisfies an evolution equation that it is convenient to write in integral form: 

^2 



^k + SjkM\v\- 



+ eij^xkxiv^j ■ V:,^ ] (16) 



where is any smooth, compactly supported test function and sij is the usual antisymmetric 
tensor. This follows from the fact that Jv = ^curlj(u) together with the equation for the 
evolution of j{v), which is obtained from ()NHGP after multiplying by and rewriting the 
result. 

Identity ()16p is central to our analysis of vortex dynamics, as in previous works [5} [6} 1151 
[21 [H] on the homogeneous case (for which of course ()16p still holds, with jy = 1). Under the 
conditions that Jv is approximately a measure of the form vr Yl\=i di^£_i(t)-, where ^i{t) are the 
vortex locations and di € {±1} their signs, one expects the left-hand side of (fTBI) to satisfy 



Jt lo^'^'' ^ dt j^^^^H^Am) - ^(vr5Zd*9'te(i))) = 7^5^diV(^(^,(^))•e^(^)• 
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Assuming that this holds, then to understand the vortex velocities ^j, it only suffices to 
understand the right-hand side of (jl6p . It turns out that it also suffices to consider test 
functions that are linear near each vortex. For such test functions, in the homogeneous 
case Vrf = 0, the integrand on the right-hand side of (jl6p is supported away from the vortex 
locations, and one is thus able to control vortex dynamics by controlling terms of the form 

■ away from the vortex cores. This argument is a key feature of all existing work on 
vortex dynamics in the homogeneous case. 

When V?7^ 7^ 0, it becomes necessary to control terms like -Vx^ across the vortex cores. 
Carrying this out, in particular relying on the approximation given by Proposition [3l is the 
main new point in our analysis. Once this is established, the whole argument is completed 
by using a Gronwall type argument on a quantity related to r* , namely T,^ + g{rl^), where 
the function gh is defined in (I19p . This demonstrates in particular that the new information 
found in Proposition [3] is strong enough to close the estimates and conclude the proof. 

Acknowledgments. This research was partially supported by the National Science and 
Engineering Research Council of Canada under operating Grant 261955, as well as by the 
projects "Around the dynamics of the Gross-Pitaevskii equation" (JC09-437086) and "Schro- 
dinger equations and applications" (ANR-12-JS01-0005-01) of the Agence Nationale de la 
Recherche. 

2 A useful lemma 

We frequently use the norm. The specific convention we use is in our definition is 

ll/^llvK-i.i(o) := sup{(;U, if) : if e Wo'°°(f]), max{||(/?||oo, ||Vv3||oo} < !}• 

In this paper, we will only use this norm on measures or more regular objects, although of 
course it is well-defined for a somewhat larger class of distributions. 

The following lemma, which we will use numerous times, is an easy special case of classical 
results (see [3] for example). 

Lemma 1. Suppose that ft is an open subset o/M", and that {ai}[^i are distinct points in 
il. Define pa ■= min{{i|aj — aj\}ijAj U {d{ai,dQ)}i U {!}}. Given any points ^'^ ^ 

and {dilLi e {±1}'. if 

' 1 
II ^(ii5(ai - 1(0) < -^Pa, (17) 

1=1 

then (after possibly relabelling the points 

I I 

i=l i=l 

In the remainder of this paper, we will always tacitly assume that under the conditions 
of the lemma, the points are labelled so that the conclusion holds. 
We give the short proof for the reader's convenience. 
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Proof. For i = I,. . . ,1, define ipi{x) := di{\pa - \x - ai|)+. Then max(||(/9i||oo, ||Vv3i||oo) = 
max(ipa, 1) = 1, for every i, so 

I I 

i=l j=l j 

Then (fT7|) imphes that {^j}^-^^ n B{ai, pa/2) is nonempty for every i. Since {B{ai, pa/2)}l^i 
are pairwise disjoint, it follows (after possibly reindexing) that {^j}j^inB{ai, Pa/2) = for 
all i. Now let ip = cpi. The functions {ipi} have disjoint support, so max(||(/9||oo, IIV93II00) = 
1, and thus || Y^\=idid{ai - Ci)\\w--^'^iQ) > {J2i=idi{6a, - = ELi " the 

other hand, if ip is any compactly supported function such that max(||'(/'||oo, || V^||oo) < 1, 
then 

I I I 

1=1 i=l i=l 

Hence || Yl\=i di6{ai - Ci)\\w-^^^n) < E!=i " ^^l- D 

3 Relating weighted and unweighted energy 

In this section, we relate the weighted and unweighted energy under some localization as- 
sumptions on the Jacobian. For a measurable subset ^4 C and v S H^{A,C) we set 

Ee,ri{v; A) := / ee,r^{v) and Ee{v] A) := Ee,i{v, A). 
J A 

Define the function g on by 



. I log x\ ■ r , 1 

X + 1 if X > Tj-!- 

otherwise. 



5(^) = <1 I+loJfe1| . (19) 



We have 



Proposition 1. Let ^} C an open set, {ai}\^i distinct points in Q,, {di}[^^ E 
and ?7 : — > M a positive Lipschitz function such that inf q 77 =: ij^in > 0. Set pa '■ = 
min{ { ^ (i(ai, flj)} i^j U {d(aj, (9r2)}j U {!}}, antf let e < exp(— ^) and v G i?"'^(r2,C) &e suc/i 
that ^ 

^a:=(^f^-7rj2vHaA <+oo. (20) 



i=l 

Assume also that 



loge| 

t — ± / 

I 

Pa 



Ta := \\Jv - vr^(ii5aJlH/~i,i(n) < (21) 
T/ien i/iere exists a constant C , depending only on I, ||V77^||oo and rjmin, such that 

< vr + C(Sa + c/(ra)) fort = 1,..., I 

ioge , , 

f22l 

E.>;ii\uUB(„..i;)) ^ 

|loge| 



8 



where R = 4max(ra, |loge| ^) < ^ and e := , and the function g is defined in ([12]). 

Proof. Let r G [r^, ^] be a number that will be fixed later. Then the balls {-6(0^, 4r)} 
are disjoint and contained in O. Let i € {1, • • • , /}, by monotonicity of the norms with 

respect to the domain, we deduce from (f2T]) that 



\\JV - 7rdi6a^\\w^l,i(^B{a„4r) < < r. 

It follows from the lower bounds estimates of Jerrard [8] or Sandier [18] that 

4r 

Es{v,B{ai,4r))>7Tlog—-Ki, (23) 



for every 5 > 0, where Ki is a universal constant. We next write 



4e2 



>/ ^'(^)[^ + ^^ tt] (24) 



> inf r/2(x)^e(t^,S(ai,4r)). 

V2:eB(ai,4r) / 

Therefore, from ()23l) with the choice (5 = e, and noting that | logr| > log |^| > log 8 > 1, we 
obtain 

^e,^(t;,5(ai,4r)) > r/2(ai)7r|log e| -i^2(r|loge| + |logr|), (25) 

where K2 depends only on HVr^^Hoo and r\min- 

On the other hand, we deduce from (1201) and (1251) that 



E,,^(7;,B(a„4r)) < S,,^(t;,f]) - E,,,(i;, S(aj, 4r)) 

i+i (26) 
< 7rr?2(ai)|loge| + S^lloge] + (/ - 1)^2 (r|loge| + |logr|) . 

Hence, going back to (I24p we obtain 

E^{v,B(a,Ar)) < 7— ^—^E,Jv,B{ai,4r)) 

(infxGB(a„4r)r(3;)j (27) 

< vr|loge| +i^3(Sa|loge| +r|loge| + |logr|), 

where depends only on ||V?/2||oo and ijmin- 

Concerning the energy outside the balls i?(ai,4r), we have from (I20p and (1250 

E,^^{v, n \ Ui5(ai, 4r)) = Es,^{v, n) Ee,r,{v, B{aj,4r)) 

i 

< Sjlogel +ZK2(r-|loge| + |logr|). 



(28) 



Hence, 



Ee{v, n \ UiB{a„ 4r)) < ^^E,^^{v, ft \ U^B{ai, 4r)) 

mfr/^ (29) 

< ir4(Sa|loge| +r|loge| + |logr|). 
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where K4 depends only on I, ||V?/^||oo and r]min- 

The function r 1— )■ r + | log r|/|log e| is minimized taking r := max(rQ, ^ , ), in which 
cas r < ^ by assumption on Td and s. The conclusions fj22|) follow with the choice C '. — 
max(il'3, 7^4). □ 

Remark 3. If we define := "^[^g^^^j^ ~ ''^ Yli=i ''f i'^i) > then (j25]) implies that 



|loge| 



for every r £ [r^, Choosing r = max(ra, ) above, we find that > —lK2g{ra)- 
In particular, = (Sq)"^ < + 2lK2g{ra)- So all our estimates remain true if we replace 
C{Ea + g{ra)) by C{ta + {2IK2 + l)g{ra))- 



4 Improved localization for Jacobians 

In this section, we prove that if the Jacobian of a function v is known to be sufficiently 
localized, then, provided the excess energy of v with respect to the points of localization is 
not to big, the localisation is actually potentially much stronger. A result in the same spirit 
was obtained by Jerrard and Spirn in [10] for the Ginzburg-Landau functional without a 
weight. Our proof here below makes a direct use of Theorem 1.1 and Theorem 1.2' in [10] by 
relating the weighted and unweighted Ginzburg-Landau energies according to Section [3j 

Proposition 2. Let Q. <ZM? be a bounded, open set, {oj}'^]^ distinct points in Q., {diY-^^ S 
{±1}, and r/ : — )• M a positive Lipschitz function such that inf^r/ =: r/min > 0. Set 
Pa = min {{i(i(aj, aj)} U {d{ai,dQ)}i U {1}}. Let £ < exp(— ^) and let v G H^{Q, C) be 
such that ^ 

^a:=(^f^-7rj2vHaA <+oo. (30) 



llogel 



Also, assume that 



I 

Ta = \\JV -TT^di5a,\\w-^^{U) < (31) 
i=l 

Then there exists Ci > 1, depending only on a lower bound for pa and rjmin o.'^d on an 
upper bound for I and ||Vr7^||oo, and for each i G {1, • • • , /} there exists a point S,i G B{ai, 2ra), 
such that 

I 

\\Jv - vr^(ii%||iy-i,i(n) < = r^{T.a,ra) = eexp(Ci(Ea + 5(ra))|loge|). (32) 

i=l 

where g is defined in Proposition [J[ 

Remark 4. Note that Lemma[J\ and (j31|) . (j32p imply that 

I 1 

Vki-Cil < -(ra+rg). (33) 



1 = 1 
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Remark 5. Since g{r) > for every r, our requirement that Ci > 1 implies that 

> e|loge|. (34) 

As mentioned, the proof of Proposition[2]relies very heavily on estimates from |10] . Follow- 
ing the proof, we discuss some small adjustments we have made in employing these estimates 
here. Also, from here upon in many places we will denote by C constants whose actual value 
may change from line to line but which could eventually be given a common value depending 
only on I, pmin, Vmin and ||V??2||oo- 

Proof. Since e < exp(— our assumptions imply that the hypotheses of Proposition [T] are 
verified. Then, since B{ai, ^) C B{ai, R)U{n\ u\^^B{ai, R)) for any R< ^, and recalling 
that g{r) > for aU r, we deduce from (I22p that 

log(^) |loge| |loge| 

for i = 1, . . . ,1, and similarly (j22p implies that 



E,(.;0\uU|^(a.,^)) ^ 

|loge| 

According to Theorem 1.2' in [TOj, it follows from (|3ip and ()35p that for every is{l, ...,/}, 
there exists some S^i S B{ai, 2ra) such that 

||Ju - 7r(ii(55J|^y-i,i(B(a^,i^)) < C eexp[C{^a + g{ra))\loge\]. (37) 

In addition. Theorem 1.1 in [10] implies that if V is any bounded, open subset of then 

\\Jv\\w-i,i^V) < C e Eg{v, V) exp ^ ^ < C e exp {Eg{v,V)) . (38) 

In particular, this and ()36p imply that 

\\J'^\\w-^-Hn\u\^,Bia,if)) < Ceexp[C(Sa + ff(ra))|loge|]. (39) 

Now for i G {1,...,/}, let Xi ^ C^{B{ai,^)) be functions such that Xi = 1 on B{ai^), 
< < 1, and llVxilloo < Cp-\ Also, let xo = 1 - E\=iXi- Then for any if G Co^(J7), 
such that ||v'||vFi'°°(n) — 1) 

« I I 

{ip,Jv -■K^diS).^) = ^{xj'^,Jv - -K^did^^) 

4=1 j=0 i=l 

i=l 

I 

< X] IIXiV'llH/i.'-C'eexp[C(Sa + 5f(?^a))|loge|) 
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where we have used (f37I) for i = 1, . . . , / and (f39]l for i = 0. Thus 



{ip,Jv -TT^^did^^) < —eexp[C{^a + g{ra))\loge\] < eexp[C7i(S<, + 5(ra))|loge|] 



for a suitable Ci, dependmg on the lower bound po for pa as well as l,r]min, ||V7y^||oo- This 



To facilitate comparison between some facts that we have used above and the precise 
statements in [10], we make the following remarks. 

First, we have used some estimates in cruder but simpler forms than they appear in |10j . 
For example, on the right-hand side of (|37|) . we have replaced an expressions of the form 
(C + Kq)"^ exp{^) from [10], where here we take Kq = C(S(j + g'(ra))|loge|, by the simpler 
expression Cexp(/Co)- We have also used the fact that Kq = C {T,a + g{^a))\^og e\ > log |loge| 
to allow us to absorb some lower-order terms from |10j . 

Second, estimates in [10] are stated in terms of a slightly different norm, ||/^||^-i.i(y) := 
sup{(//, (/>):(/)£ C^{V), \\V4>\\oo < 1}- This does not cause any problems for us, since clearly 

ll/^llvK-i.i(v^) - llA'lliy-i>i(v^)- 

Finally, the estimate corresponding to (|38p in [10| is a special case of a more general 

result, and as stated there requires the additional assumption that "^log^^ < tt. However, 
since || Jf < ||'/'w||Li(y) ^ '^Eg{v;V), it is clear that (j38|) is still true if "^ipg^^ > tt. 

Remark 6. IfQ is simply connected, then we can alternately argue by citing a result from 111^ 
to obtain an estimate of the form (I32p with Ci independent of pa, at the rather small expense 
of having to replace ^ on the right-hand side of (jSip by some smaller quantity depending 
on I as well as pa- This is in principle useful if one wants to consider large numbers of 
vortices. The relevant result (Theorem 3) from 11 If is proved using facts from UOf . as in the 
proof above, but combining estimates on the balls and away from the balls in a more careful 
way, to avoid introducing the factors of p~^ that arise from the cutoff functions that we have 
employed here. 

The proof of Theorem 3 from ]llf can surely be adapted to yield a similar result without 
the assumption that Q, be simply connected, but since the proof is slightly complicated, we 
prefer not to tinker with it here. 

5 Across the core approximation by reference field 

In this section we prove 

Proposition 3. Let 17 C be an open set, {^lY^^^ distinct points in 0, {(ii}'^^ G 
and : — 7- M a positive Lipschitz function such that infj^ rj =: r/min > 0. Set p^ = 




mm {{^d{^i,£,j)}i-tj U {d{^i,dn)}i U {1}}. Let e < exp(-|-) and let v G H^{n,C) be such 



implies ([32]) . 



□ 




that 




(40) 
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and 



\Jv - 7r^di5gJ|^y-i.i(n) <r^ = eexp(K|loge|) = e 



l-K 



(41) 



1=1 



for some K < ^. 

Define = j*i{(,i},r^) in Q by 



di i^'f'-*^. iV2 ifxe Biii, TT^) 

,/3;eo\uUi?te,^) 



where {yi,y2)^ := i-y2,yi)- Then 



1 





j{v) 


[ 


1 1 J* 


In 


\v\ 



< {Ci:i: + K)\\oge\+C\og\\oge\ 



and 



(42) 
(43) 



II V X {j{v) - J*) 11,^-1.1 < C'3?'e- 
where the constant C depends only on I, ||V77^||oo o-nd rjmin- 

Since -ft' < |, the assumption that e < exp(— 8//3g) imphes that < < In 

particular, the balls -B(^i, |loge|~^), i = 1, . . . , / are pairwise disjoint and contained in il. 

Proof. We will use more than once the fact that 



|VvP = |V|t;| P + 



j{v)\ 



(44) 



Step 1: verification of ()43p . A direct calculation, using the definition of j*, shows that 
for any smooth (p, 



i=l ^5 JB{£,,,rf^) 



Thus ||V X j^, — 27r ^ (ij(5^J|i4/-i,i(Q) < Cr^. Recalling that Jv = x j{v), we deduce 
from this estimate and our assumption (j4ip . 
It remains to prove (|42p . 

Step 2: decomposing the energy. Note that our assumptions (j40p . (j4ip about the 
points are exactly the same as the hypotheses ([20]) . (f2T|) about the points {a,}^^]^ 

in Proposition [H except that here we impose an additional smallness condition on r^. Thus 
estimates from Proposition [T] are all available here. In particular, recalling (I28p with the 
choice r = max(rg, i ' 



li^l^, we see that 



E,^r,{v,^\^iB{i,,- -)) < C7(S^|loge| +log|loge|). 

|loge| 

In view of (j44p . and noting that is supported in Uji?(^j, ) to prove (|42p it therefore 
suffices to show that 

2 



log e 



3* 



< ((:7S5 + K)|loge|+(:71og|loge| (45) 
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for i = 1, . . . ,1. Toward this end, note that 

, log I log e 



|loge| 



7rry^(Ci) + S5 + C- 



|loge| 



> Ee,r,{v,B{C,,- -)) 

\loge\ 









1 1 




\v\ 



1 



J* j 'J* 



Using the explicit form of j* and of , 

1 



?? > - mm T] 



> (??'te)-C||Vr/'||oo|loge|-^)7rlog 



|log£| ^ 



By combining the previous two inequalities and rearranging, we see that to prove (j45p . it 
suffices to check that 



<C(S5|loge|+log|loge|). 



(46) 



Step 3: proof of ()46]l . 

First note that 



2/j(^) • N • 

V (tt - J*) - J* 



+ V\^i) [ 
Je 



j* (1 - l^^l 



- J*) - J* • 



We estimate the three terms on the right-hand side in turn. First, 



< C 



Ee,rj{v) 

|loge| 



+ (1- J^+|loge|-i: 



<C(S5 + log|loge|), 

where we have used the fact that |log e|"^£'e^,,(u) < C(S^ + /7r||?7^||oo) < C(S^ + log |log e|). 
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Next 

„2 



•j; (i-IH)l <c|b;||oo / ^|vH' + :^(bp-i)' 



< C(S5 + log|loge|), 



using the lower bound for and arguing as above. 
To estimate the final term, note that = V^/i, for 



(di[^{\x-Ci\^-l)+\ogir^\loge\)] if\x-^,\<r^ 



h{x) := < 



di \og{\x - Ci\ |loge|) 




if < \x — < |log e| ^ 
if X 0Ui5(e»,|loge|-i). 



Thus, we can integrate by parts to find that 



{j{v) -j*)-j*\ = \ hV X {j{v) - j;)| 

< max(||/i||oo, ||V/i||oo) \\j{v) - j*||vi/-i.i 
<C, 

after using (j43p and noting that ||V/i||oo = ||j||oo = f'^^- This completes the proof. □ 
It follows from the definitions ([30]) and (jlO]) of and , together with (p3|) , that 

< + C<7(ra) 

for C depending only on ||V?7^||oo- Combining this with Propositions [2] and [3l we immediately 
obtain 

Corollary 1. Under the assumptions of Proposition\^ 

2 





j{v) 


[ 


1 1 J* 


Jn 


\v\ 



< C(S„ + 5(r,))|loge|, 



(47) 



where = j*{{£,i},r^), the points are given by Proposition and the constant C 

depends only on I, pa, ||Vr/^||oo and rymin- 



6 Small time upper bound on the speed of vortices 

Let Ci be the constant given by Proposition [2] corresponding to the lower bound Pmin 

(as 

defined in (161)) for pa-. Let also e < exp( —)■ Then the conclusions of Proposition [2l 

pmin 

applied to with this choice of constants, are available to us for all < t < T^oi- Since the 
conclusions of Proposition [2] remain true if we increase Ci, we may assume that 



1 

— < 



(48) 
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which we do in the sequel. We define the stopping time 



Tioc = sup{t < T^oi ; + g{rl) < V < s < t}. 

Since the function g satisfies g{r) > r on M^, for t < Tioc we have r* < 2^ < ^"g" . In 
particular, we may apply Proposition [2] to w*, {ai(t)}'^^ and {dj}^^]^, which yields points 
{ii{t)} such that 

llJ^t _ ^^d,%(,)||H/-i,i(n) < 4 = ^c(S*,r*) = £exp(C7i(S* +5(r*))|loge|), (49) 

i=l 

wher^ 

Since t 1— )■ is continuous in H^{Q,), it is clear that t 1— )• Ju* is continuous as a function 
from M into and hence we can choose {Ciii)} to be piecewise constant, and in 

particular measurable, as functions of t. Since E^^^] is preserved by the flow for v and rj^{ai) 
is preserved by the flow for the a^'s, we have 

Note in particular that r| < y/e for t < Tioc- 

Proposition 4. There exist positive constants ro,eo ond C, depending only on I, pmin, Vmin 
and ||Vr7^||oo; such that Eq < exp(— — ^) and if < e < Eq and 

Pmin 



4Ci 

for some t < Tioc, then Tioc ^ t + tq o-n-d 

\\Jv'-Jv'\\w-i,nn)<C{\t-s\+rl), (50) 
r| < r| + C|log e\e^/^{\s - t\ + r^) , (51) 
{a,{s),^,is)} C B{a,it), ^), i = l,...,l (52) 

for every t < s < t + tq. 

Proof. For the ease of notation in the present proof, || • || is understood to mean W~^'^{Q) 
while I • I denotes the Euclidean norm on M^. 

Step 1. Let t < s < mm{Tioc,t + tq}, for tq to be fixed below. We first use the fact 
that Jv^, Jv^ are well-approximated by sums of point masses to show that \\Jv^ — Jv^\\ can 



^Proposition [2] actually uses a version of surplus energy for which the weighted energy Ee,n is restricted to 
Q. Since the energy density and the weight are non-negative, our definition of surplus E* here, integrating on 
the whole R^, yields a larger number, and is therefore compatible with the claim of the proposition. 
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be estimated by computing (Jv^ — Jv*,ip) for a specific test function ip with certain good 
properties (in particular, bounds on second derivatives of ip). Toward this end, note that 

I I I 

\\Jv' - Jv^W < \\Jv' - TT'^di6^^^s)\\ + 11-^^^* - ^X^'^^^W" + 

'=\ (53) 

<rl+rl + 7rj2Us)-m\- 

1=1 

We now fix tq, depending only on ||V7/^||oo, Vmin and pmin, such that if t < s < t + tq, we 
have \ai{s) — ai{t)\ < ^^^^ for all i G {I,-"" J}- By Proposition O the choice of Tioc, and 
Lemma dl for every r < T/oc we have |aj(r) — Ci(r)| < 2r2^ < £s^^ gy i^j^e triangle inequality, 

4 



it follows that ^i(s) e B{ai{t), for ah t < s < min(t + tq, T^oc) and i e {1, • • • , Let 



( \ ^ .J {x-ai{t)) ■ iCi{s)-Ci{t)) /, .,^|^ frA\ 



i=l 



where x e C'^(M+, [0, 1]) is such that x = 1 on [0,pr«m/4], x :se on [^^f^.+oo). By con- 
struction and the definition of pmin, we have ip £ I?(0) and it follows that 

I I 
^^Ms) -£.iit)\ = {TT^di{5^^(^s) - Sut))^'P) 

i=l 1=1 

< (r^ + r|)||v9||^i.oo + {Jv' - Jv\ip). 
Combining this with (I53p . we conclude that 

\\jys _ j^;t|| < C(r| + r*) + {Jv' - Jv\^), with \\ip\\w2,o. < Cp^^^-^. (55) 

Step 2. We now deduce from (j55]) . together with (fT6|) . the fact that S'' < and 
conservation of energy, that 

II J-u"* - Jw*|| < (r^ + r|)||(/p||y(/i,cx) + (s - sup ||-^ Jf^||vy-2,i(n)ll¥'llvy2,°o 

rG[t,s] (56) 

<C(r| + r| + |t-s|), 

for t < s < min(t + to,Tioc), where C depends only on /, pmin, Vmin and ||V7?^||oo- 

Step 3. It remains to estimate r| and to show that t + tq < Ti^c- For that purpose, since 
r| = r^(S'^,r^), we first use (1561) to compute 

I 

i=l 

' ' (57) 

< llJt;' - Ji;*|| + \\Jv^ -7r22diSa,{t)\\ + h22^i(^^d^)~^^dt))\\ 

1=1 1=1 

<ri + C{\t-s\+rl + rl)). 
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Next, since s < Tioc, 



r 



■l = eexp{Ci[E^ + g{ri)]\loge 



<r\ + Ci\\oge\e-2\\g'\\o^ {ri-r\)+ (53) 
1 



provided we assume, and this is again no loss of generality, that C|logeo|eo ^ 2U 
same constant C as in (j57p . Combining (|57p with (j58p we obtain 



<-r* <C(|t-s|+r*). (59) 
Going back to ([58|) . this yields the desired estimate of r| : 

r| < r\ + C\\oge\e^''^{\s-t\+r\). 

Then going back to ([56]) , 



II Jz;' - J^;*|| < C{\t-s\ + r\) 

ioi t < s < min(t + tq, Tioc)- Finally, by assumption we have + ) < l/(4Ci) so that by 
([591) and the fact that g' <l, 

S° + 5(r-^) < l/(4Ci) + C7(|t - s| + r*) < l/(4Ci) + C7(to + e^) < l/(3Ci), 

1 

provided we assume, and this is no loss of generality, that C{tq + £q) < l/(12Ci). It follows 
that min(t + tq,Tioc) = t + tq, and the proof is complete. □ 

7 Control of the discrepancy 

In this section, we prove a discrete differential inequality for the quantity r* . More precisely, 
we will prove 

Proposition 5. There exist positive constants eq and Co, depending only on I, pmin, Vmin 
and ||Vr?^||oo, such that eq < exp(— — ^) and ifO < e < Eq and 

Pmin 

S°+5(r*)<^ (60) 



for some t < Tioc, then 

<Co(SO+5(r*)) 



T-t 



where T = t -\ ^ < Tioc- 



This is the main estimate in the proof of Theorem [2j 
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Proof. We first require the constant eo to be smaller than the one appearing in the statement 
of Proposition SI As in the proof of Proposition HI we will write simply || • || to denote the 
^^'^{Q) norm. Note that the condition ()60|) states exactly that 

4 < e^/^ (61) 

and then the definition of T and ()5ip yield 

r| < 2r| for ah s £ [t, T] (62) 

if C is large enough and eq small enough, which we henceforth take to be the case. Moreover, 
from ([59|) . we see that < rli^ + C{T — t + r*^) for all s G [t, T], and then the choice of T and 
the definition of g imply that 

g{r'J<2giri) for ah s S [t, T]. (63) 

1. First note that 

I I 

i=l i=l 

I 

<7^Y,mT)-ai{T)\-Ut)-a^it)\) + rj + r\ 

i=l 
I 

< vr ^ i/i • {UT) - m + a^{t) - ai{T)) + rj + r* (64) 

i=l 

for fi = (unless ^i{T) — ai{T) = 0, in which case Vi can be any unit vector). We 

now define 

^{x) = '^dii^i ■ - ai{t))xi\x - ai{t)\) 

i 

for X S C°°(1R"'", [0, 1]) such that x=^on [0, ^Pmin] and x = on {pmin, oo). It follows from 
(l52l) that (since d,-^ = 1 for ah i) 



vr ^ I/, • (^i(r) - ^i{t) + a,{t) - ai{T)) = i^Y^di [p{UT)) - ^{Ut)) - ^{ai{T)) + ^{ai{t)) 

i=l i=l 

so that (|64p and the definition of imply that 

I 

rl-ri<{^,Jv^ -Jv') - TTY,di[p{ai{T)) - ^{ai{t))\ + C(r[ + r*). (65) 

i=l 

2. The remainder of the proof is devoted to an estimate of (93, Jv^ — Jv^). First, using 

fT Q 

{ip,Jv^-Jv*)= / —((p,Jv'')ds 
Jt ds 

I f'^ f f (beP - 1)2 \ 

J J {^Ijy^xtVl^ ' 2 +(^ljVxkX,Ve,x, ■Ve,x^j (66) 



|loge| 



+ / / ^'i'^^iT;2 



t Jn 



rf |loge| 
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We immediately see from ()47p that 

1 



|loge| 



2 i\v\'-l? 



<C(S° + 5«)) 



for every s G [t,T]. Moreover, it follows from ^2i) that 4|log e|"^) C B{ai{t), \p, 

if eo is small enough, and the definition of 93 implies that i^XiXj = in u'^^-B(aj(t), ^/3- 
so (j'22|) implies that 

^ • v,,x,\ < C^^^Ee{v- n \ UB{^u 4|loge|-^)) 

<C(S0 + 5(r|))<C(S0 + 5(r^)) 

for every s G [t,T]. 

3. We now decompose the remaining term in (I66p . For every s G [^j^]; l^lH 
j^=({.^j(s)}, be the approximation to j{v^) obtained in Proposition [3l Note that 



(67) 

min) 
'mm) ) 



V V 



where for example j{v)j = {iv,dx v) denotes the jth component of j{v). Thus, adding and 

subtracting in various places, and writing •i/'^-fc as an abbreviation for €ij(pxi^T, we have 
for every s G [t,T], 



'^e,Xj ' '^£,Xf^ 



|loge| 



|loge| 



+ / ^jfc 

+ / Ipjk 



1 



|loge| 
1 

|loge| 



jl if*)k + 



j* I U*)j 



(69) 



I I * 



We immediately dispense with the easiest terms by using (j47p to see that 



/ V'ifci r 



3 ^ 



<C{TP + g{ri)) (70) 



for every s G [t, T]. 

4. We next consider the first term on the right-hand side of ()69p . which is the term that 
yields the dominant contribution. Since jl is supported in Uj-B(^j(s), |loge|~^), clearly 



U*)j{j*)k 

|loge| 



E 



if^jk 



{jt)jU*)k 



^yB(C,(s),|iog£l-i) |loge 



^Note that the regularization scale r| is fixed for s G [t,T] 
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For each z = 1, . . . , Z, if x G B{S,i{s), |loge| then \x — ai{s)\ < |loge| ^ +r^ + r|, by 
so for every s £ [t,T], 



U*)j{j*)k 



lb* 111 

|loge| 



'B{6(s),|loge|-i) ^ ^ I^Ogel 

< C7(|loge|-i+< + r|), 
using the expHcit form of j*, which (together with the definition (I49p of r|) also impUes that 

ys(5.(s),|iog£|-i) |loge| |loge| r| 4 

= .5,.(l-C,(S°+,(rD)) + O(i^lll^). 

Combining the above computations and recalHng that g{r) > max(r, i2i-H2I£i) for all r and 
that g{r^a) ^ ^| ^ ^ioo we conclude that 

Ja |loge| ^ 

= +0(C^i(S° + <7(r:)). (71) 

In the last line we have used the definition -i/'fefc = ^Wxidxt.i^ogrf') = V^?- V^(log 77^) together 
with the ordinary differential equation ([5]) satisfied by the points ai{-). 

5. Combining and imd recalling dS^]), ([Ml), we find 

that 

rl-ri<C{T-t){i:^ + g{ri)) + Crl + ^ L^^[^-jt) {f.hdxds (72) 



t Jn 



i'jk f i{v) 



t Jn |loge| \\v\ / fc 

We now begin to control the integrals on the right-hand side above. We will consider only 
the first, since the estimate of the second is identical. First, 



3l] {iDjdxds. 




t Jn 



\loge\ V \v\ V 7 * Jn |loge 



a Jt 



|loge| \ \v\ 



3 



We claim that 

^^^^^*"^*^'^-(^"^*)- '^''^'^ <C(T-t)(SO + 5(r*)). (74) 
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Using the Cauchy-Schwarz inequality, dH]), and (163]) . we see that it suffices to prove that 

/ |J* - Jtl' dx < (S" + <7(r* ))|log £\ for every s G [t, T]. 

Jn 



Toward this end, we fix some such s, and we introduce the notation 

I 



1 ' 



Our choice of T and dSO]) imply that a < C-^. Writing Bi := B{ai{t), ^^), we deduce 
from (|52p and the support properties of that 



/ \jt-f.\'dx = Y,I \jl-jl?dx. 



For each i, B{^i,a) C B{^i{t),2a)riB{^i{s),2a), so by an explicit computation, and recalling 
and the definition of r|, we find that 

/ \jl-ffdx<2[ \ji\'dx + 2 [ |j-|2dx<21og(^) + 21og(^) + C 

< Clog(^) 

<C(SO + 5(r*))|loge|. 

Next, on B{^i, ^|^^^^| ), the definitions imply that both and f are nonzero, and in fact 



\x-m\'\x-m\'' 

Since \(,i{T) — fi| < | for r = t, s, it follows that 

Ix-^iT 2|loge| 

and hence that 

/ \jt-jtfdx <C. 

Finally, 

/ |j* -J*P dx <2 [ Iff dx + 2 [ Iff dx < C. 

•^^ABte,^^) -^sABteW,^) JBABte(^),^) 

We deduce (fTlI) by combining the previous inequalities. 
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6. We now consider the first term on the right-hand side of ()73p . Clearly 

-iji)k(4^-jiv)] dsdx 



3 i=l-"^i 



|loge| 



i=l 



^3k 

|loge| 

(75) 



By elementary estimates, 



\j{v)\ 



\v\ \v\ I Ze 

and from the definitions, and recalling (j34p . we see that ||i*||oo ^ {'<'\)~^ ^ (^llog^l)""'^- Thus 
for every i, 




T 



Bi Jt 



|loge| 



{jl)k ■ dsdx 



< C\\jl\\ooe{T-t) 



|loge| 
+ C 

<C{T-t){TP+g{rl)), 



|loge| 



since 



(77) 



i=l 

3 



7. Now fix some i G {1, . . . , Z} and let G C^iB{ai{t), jPmin)) be a function such that 
X* = 1 on Then for every s G [t,T], 

UJiv) - = V/^ + 4v V in B{ai{t), ^p^mn) (78) 

for and (7*, real-valued functions on B{ai(t), ^Pmin), solving 



(79) 



V • iv'vn = V • {rv\jiv) - jt)) in Bia,{t), ^P-™), 

3 

z^-V/'' = on dB{ai{t),-p^in), 

and 

3 

/ = on dB{ai{t), -pmin)- 

Indeed, if we let be a solution of ([79]) . then t]'^ {xiU i^) — jt) ~ ^f^) is divergence-free and 
hence can be written as V^g^ on B{ai{t), \pmin)i so that ([78]) holds. Then it follows from 
(j79p that satisfies the equation in (|80p . and that the boundary condition is satisfied after 
adding a constant to . 



(80) 
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Thus 




B, Jt 



|loge| 



{jl)k - fDidsdx 



I, r(j:)fc(VF+^),dx (81) 
|ioge| r/^ 



for 



F{x) 



G{x) 



g%x) ds. 



We write F = Fi + ■ ■ ■ + F^, where 



for 



V-{-rfVFm) = Am in B{ai{t),-pmin), • VF^ = in dB{ai{t), -p„ 



Ai = x' I V • iv'jiv)) ds 

rT 



A2 
^3 
^4 



^ ri^Vt ■ - 1) ds 



\v\ 



3l) ds. 



Using the continuity equation p5|) — this is a key point in our argument — and (j77|) . we 
note that 



L2 



\\A,\\L2 = \loge\ x'v'{\v\'-i)\^ 

<Ce|loge|3(sO + g(r*)) 
since g{r) > ^+*°s |i°g^l ^^j. Next, the definition imphes that V • j| = for every s and 

log £ 

1— - 

that IIj^IIlp < C'pllogej p for every p < 2, so 

P2IILP < WXi IIlco (T - t) sup ||V(7y2) . < C{T - t)\\oge\^"l for p<2. 

se[t,T] 

Very much as in (|76p . we can check that 

Pallii < C(T-t) sup \\^-^{\v\-l)U.<C{T-t)e\\ogE\\TP + g{ri)), 



and it follows from ()47p and (j63p that 



P4||l2 < C(T-t)(|loge|(SO+5(r*))) 



1/2 
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Clearly, for any . . . , ^4 G [1, oo] 

4 

' II A II IIT7 ), ,, , 

mllg; 

where — — h = 1. Using elliptic estimates and Sobolev embeddinff theorems, and taking 
qi = |, 

^^lb;|||l|VFi||4<^^||j;|||Pi|b < Ce|loge|i(SO+5(r*))<C(r-t)(SO+5(r*)). 



The last inequality follows from the choice of T and (|34p . which imply in particular that 
T — t > e|logep. Similarly, taking q4 = 4/3, 

^||j;|UI|VF4||4 < ^^P4||2 < C^^^(SO + f7(r*))V2 < C{T-t){^^ + g{ri)), 
|loge| ^ |loge|2 llog^l 

since llogel""*^ < girl). For any g2 G taking p2 < 2 such that P2 = q2, so that 

= ^ — we find from our estimate of Ao that 

pr^llj^ll^JlVFsll,, <-^|b;||,,P2|U <C(r-t)|loger2<C7(r-t)g(r*). 
|loge| |ioge| 

And, recalling by Stampacchia's estimate that for any p G [1,2) there exists Cp such that 
llVFallp < CpllAalli, we compute, choosing qs = 3 for concreteness, 

1 „ 



IIVF3IU < l|j;||3P3l|i <C(r-t)(r*)-le|loge|(S0 + <7(r*)) 



|loge| 

<(r-i)(e|loge|)2/3($]0 + <7(r*)) 

again using the fact that > e|loge| for all t, see (j34p . Combining the above, we find that 
for every i G {1, . . . , ^} and < e < eo with eo sufficiently small, 

[ ^Xji)k-{yF),dx<C{T-t){T.^+g{ri)). (82) 
Jb, I log e I 



8. Next, 



Jb, |loge| Jb, r|loge| 



for Gm solving 



VG 3 3 

-V • ( — ^) = A'^ in B{ai{t),-pmin), g = on dB{ai{t), -pr, 



^2 



with 



A[:= l%'Vxijiv)-f,) ds, 
A'2--= rV^x'-jiv)il-^)ds, 

Jt M 

^s--=[v^x'-{^-jt) ds. 
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The terms containing G2 and G3 are estimated exactly as the terms containing F3 and F4 in 
Step 7 above, leading to 

/ -4r^\{ji)ky^{G2 + G,),dx < C{T-t){TP + g{ri)). 
JBi V \^oge\ 

For the remaining term, we invoke the interpolation inequality 

<cp;r^_MP'illi7' (83) 

for p G (1, 2) and 9 such that | = f + (see e.g. [20] Theorem 2.4.1 combined with Sobolev 
embedding theorem). To estimate the norm, we fix ^ G C^(f]), and we compute 

(C,^'i)= / (x'C,Vx(i(^)-j-))< / ||x'Clliyi,-l|Vx(i(^;)-i:)||t^-Md5 
Jt Jt 

<C{T-t)r\\\C\\w^,^ 

using (gSD and (I62D. Thus 

\\A\\w-^.^<C{T-t)r\. (84) 

Also, for every s G [i-,T], 

II V X {j{v) - j^L^ < PMl^ + II V X fjLl < CEe,r,{v) + 27tI. 

Estimating E^^^] as usual by C|loge|(S'' + g{rl^)), integrating the last inequality from t to T, 
and combining it with (I84p and (I83p . we obtain 



m\^-u.<C{T-t){rlf (C7|loge|(SO + 5(r* )))'"'. 
Then using Holder's inequality and (again) the fact that ||i,"||p/ < C7(r|)i^ for p' > 2, 

L Wk^'*^''^^^^'^''" - (C|loge|(S° + ,(r*)))^-^ 

<C(T-O|loge|-^ {J:' + g{ri)f-' 
<C{T-t) (sO + <^(r*)), 

since it turns out that 0+^ — 1 = 0, and noting that llogel^"*^ < (^(r*) for all t. Assembling 
these estimates, we find that 

/ ^(jt).(V^G), dx < C{T - t) (SO + g{ri)) . 

Now by combining this with ([72]), ([73]), ([74}, ([76]), (l82]), we finally obtain 

r* <C7(r-t)(sO + <7(r*)). (85) 

□ 



T 

V — r 
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8 Proof of Theorem [2] 

Our main result is a straightforward corollary of the discrepancy estimate proved in the 
previous section. 

Proof of Theorem [H Let Y denote the solution of the ordinary differential equation 

Y{t) = Coh' + g{Y{t))), y(0) = rO, 



where g is the function defined in (fT9|) . and let {^nlJ^o ^ discrete approximation to Y{-) 
obtained via an Euler approximation implicit in the statement of Proposition [5l Thus, we 
define 

Yo = rl y„+i = y„ + {tn+1 - tn) Co (SO + g{Yn)) , 

''n+1 • — i^n ~r 

e 

where 

rg" := r^{YP, = eexp(Co(SO + g{Y^)))\\oge\). 

Since the function f{Y) := Cq (T,^ + g{Y)^ is convex, a forward Euler approximation to the 
solution of the equation Y' = f{Y) is always less than or equal to the actual solution, and it 
follows that Yn < Y(tn) for all t. Then repeated application of Proposition [5] shows that 

ri" <Yn< Y{tn) for every n such that t„ < T^oi and S° + g{Yn) < 

Given an arbitrary t G {0,Tcoi] such that S*^ + g(Y{t)) < there exists some n such that 
t G [tn,tn+i] and r*" < Y(tn)- Then by Proposition [H see in particular ([57|) . as well as ([62]) . 

r* < r*" + - tn) + r*") < y(t) + Ce^/^^ (86) 



since the bound T,^ + g{Y{tn)) < 4^ guarantees that r|" < e'^/^ and hence that tn+i — tn < 
e^/^. It remains to bound the function Y from above. For that purpose, we notice that since 
giy) ^ y + log |loge|/|loge| for every y > 0, we have Y{t) < Y{t) where Y is the solution of 
the ordinary differential equation 

, log|loge| ~ . ~ _ 



Y(t) = Co (so + ^' ^ ' + y (t)) , y (0) 



|loge| 

The solution of the latter is explicitly given by 

m=r°+(S° + r° + i^^i^)(e^o*-l), 

and the conclusion therefore follows from (j86p . increasing the value of Cq to the value of C 
in if necessary. □ 
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9 Some properties of the ground state 

In this section we briefly recall some facts about minimizers of the functional^! 

£,y{u) = f\^ + ± (v{x)\u\' + \\uA dx (87) 



in the space 

Urn ■■= {u £ H^iR^; C) : / V\u\^ < oo, / \u\'^ = m} (88) 

where V : — )• [0, oo) is a smooth function such that V{x) — t- oo as |x| — t- oo, and m > is 
a parameter. 

For every positive e, m, the existence of a function r]s^m '• ~^ (0, oo) minimizing £i,y in 
Tim is standard, and follows easily from the growth of V (which implies that the constraint 
is preserved for weak limits of sequences with equi-bounded energy) together with the strong 
maximum principle and the fact that i^e^v'd^l) < £e,viu) for all u. 

In the introduction, we already introduced the unique number Aq such that 

(Ao — V)^dx = m, 

T 

and we have denoted by pj.p := (Aq — V^)"^ the Thomas-Fermi profile associated to V and m. 
We also note w := (Aq — V)~. We will prove 

Proposition 6. Let r] = r]^_m G "^m be a positive minimizer of £s,V in- T~Lm,- 
Then 

Wrl" - PTF\\Lm^)<Ce^'\ (89) 



Moreover, for any K CC I^tf := {x G : Ptf{x) > 0}, there exists a constant C = 
C{m,V, K) such that 

\\r]^ - PtA\l^{K) < Ce^/\ ||Vr?2||^^(^) < C. (90) 

This is quite standard, and is proved for particular potentials y in j7] for example. We 
include a complete proof, since the references we know all impose slightly more restrictive 
conditions than we consider here (for example, symmetry conditions, or the assumption that 
Ao is a regular value of V). 

Proof. It suffices to prove the result for e < £o, for some eo > 0. 
1. First, as is standard, for u G Tim we rewrite 



/ IVtil 1 /I i2 \2 1 I l2l r 1 /^ ^ 1 / 2 

£e,v{u)= I -t;— + - Ptf) +^w\u\ dx + Ta^-^o^-J / Ptf 



=: £e,pj,piu) + Ci{e,m). 



e^' " 2 4 



Thus, it is clear that a function minimizes £sy in T-Lm if and only if it minimizes £e,pTF 



Note that we make no restriction on the dimension A'' here. 
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2. Next we claim that 

inf£:,,p^^ <Ce-2/3. (91) 
Note that this immediately imphes (j89l) . We verify (I9ip by choosing '■= Csfe{\fPj,p)^ where 



feis) 



£ °'s^ if s < e*^ 
s if s > e*^ 



where is chosen so that Ue G T~im- Then straightforward estimates very much hke those in 
[7], for example, show that £e,r,iUe) < C(e-" + e^""^), and ([HI]) follows by taking a = 2/3. 
(This crude estimate has the advantage of holding for every m > 0, so that we do not require 
Ao to be a regular value of V. If Aq is a regular value, then a variant of the same construction 
shows that inf-^^ £e,pTF — C'|loge|.) 

3. Since V — Xq = {V — Ao)"*" — {V — Ao)~ = w — Ptf, we may write the variational 
equation satisfied by rj in the form 



-Arj + ^iv^ - Ptf + w)r] = ^(A. - Ao)r?, 



where ^A^ is a Lagrange multiplier. Multiplying by rj and integrating, and using the fact 
that T] G Hm, we find that 



m 

It follows that 



{Xe - Ao) = / |V7?|2 + 1 [wT]'^ + {if - p^pf + {rf - Ptp)p^ 



^{Xe - Ao) < 'i£e,pTFir]) + -^\\PTF\\L2{m)\\v'^ - Ptf ||l2(ik]V) < Ce (92) 



£2 



by and ([Ml). 

4. Now let pTF,e ■■= (Ae - V)+. It follows from ^ that 



PTF,e — Ptf\\l°°(m.'^) — |Ae — AoI < Ce^'^^, (93) 



so that K CC := {x G : Ptf.e > 0} if e > is sufficiently small, which we henceforth 
take to be the case. Note also that 

- AT] + \{ri'^ - PTF,e)r] = mQe- (94) 

Now fix some r < ^ dist{K, dQ,TF)- In view of (j93p . and since V is C^, there exists a, A; > 
and £0 > such that 

Ptf.e > and |Ay/py^^| < k whenever < e < eo- (95) 

For any x £ K and b £ (0, a), define 

G,fe(y) = Ciy) = - 1)' 
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in B{x,r). Then for b G (0, f), 

1 3a^ 

-AC + ^(C'-/5tf,)C < -AC-^C < in B{x,r) (96) 

whenever e is sufficiently small. It follows that r] > (^x,b in B{x,r) for every b € (0, |), as 
otherwise we could find some 6o S (0, |) such that min5(2.^r)(r/ — Cx,bo) = 0- Since r? > 0, the 
minimum would have to be attained in the interior of B{x,r), and this is impossible in view 
of (HH) and ([96]). 
It follows that 

9a 

v{y)> — =:amB{x,r/2). (97) 

Note also that ||??||Loo(jjiV) < ||^PTF,e||L°°(R^)i since otherwise fj := min(r7, ||i//JTF,e||oo) would 
satisfy S^^prpp {fj) < £s,pTF('n)i contradicting the minimality of 7]. 
5. Now write 6 := r] — yJPTF,e- Then 

1 i3 a2 

-A0 + ae(x)^ = A^/p;^ for a£(x) = ^(0 + 2^/^77;) (^ + Vp^) > ^ 

in B{x^r/2), and \6\ < ^W^/p^^ J\iaa(^N-j on B{x,r/2). Now for y G B{x,r/2) define 

k 



I 1 2 2 

a |y — x| r , 
Ve^ 2 8"' 



where A; is the bound for || Ay^/V/TTlloo found in ([95]) . Then 9 > ^ on dB{x,r /2), and there 
exists Eq > such that 

{-A + as)@ > k > {-A + ae)e in B{x,r/2), ifO<e<eo- 

It follows that e > 61 in B{x,r/2), and similarly -G > -9 in B{x,r/2). Thus 

\r] - ^/p^\ < Ce^ onB{x,r/4). (98) 

6. Returning to we see that 

-A?/ + 6e7/ = ini?(x,r/4), for 6^ = ^(r?^ - pe), 

and ([98l) implies that ||?)£||Loo(_B(x,r/4) < independent of e G (0,eo) and x £ K. Since 
we already know that ||?7||l°o(r'V) ^ Ci we conclude from standard elliptic regularity that 
\\'^'n\\L°^{B{x,r/8)) ^ C. Also, it foUows from (f93|) and ([98]) that ||ry^ — /3||Loo{_ft') < Ce^/^, so we 
have proved ([90]) . □ 



10 Proof of Theorem [T] 

In view of Theorem [1] is a direct consequence of Theorem [2] combined with Proposition [5] 
and the continuity of the solution of an initial value problem with respect to the nonlinearity. 
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